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ABSTRACT 

We report on the detection with the MAGIC telescopes of very high energy 7-rays from IC 310, a 
head-tail radio galaxy in the Perseus galaxy cluster, observed during the interval November 2008 to 
February 2010. The Fermi satellite has also detected this galaxy. The source is detected by MAGIC 
at a high statistical significance of 7.6 a in 20.6 hr of stereo data. The observed spectral energy 
distribution is flat with a differential spectral index of —2.00 ± 0.14. The mean flux above 300 GeV, 
between October 2009 and February 2010, (3.1±0.5) x 10~ 12 cm" 2 s" 1 , corresponds to (2.5±0.4)% of 
Crab Nebula units. Only an upper limit, of 1.9% of Crab Nebula units above 300 GeV, was obtained 
with the 2008 data. This, together with strong hints (> 3er) of flares in the middle of October and 
November 2009, implies that the emission is variable. The MAGIC results favour a scenario with the 
very high energy emission originating from the inner jet close to the central engine. More complicated 
models than a simple one-zone SSC scenario, e.g. multi-zone SSC, external Compton or hadronic, 
may be required to explain the very flat spectrum and its extension over more than three orders of 
magnitude in energy. 

Subject headings: gamma rays: galaxies — galaxies: active — galaxies: individual (IC 310) 
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1. INTRODUCTION 

Most of the presently known extragalactic very high 
energy (VHE, above 300 GeV) 7-ray emitters (~40) 
are blazars. So far o nly two radio galaxies, M 87 
(lAharonian et all I2003D and Cen A (lAharonian et al 



200S ) , and two s tarburst galaxies, NG C 253 (|Acero et al 



20091 ) and M82 (jKarlsson et alj|2009f ). have been clearly 



identified in this energy range. 

IC 310 (redshift z=0.019) is a head-tail radio galaxy 
located in the Perseus cluster at 0.6° (corresponding to 
~ IMpc) from the cluster's central galaxy, NGC 1275. 
Head-tail radio galaxies display a radio morphology con- 
sisting of a bright head, close to the optical galaxy and a 
fainter, elongated tail. In the standard explanation, the 
jets are bent towards one direction creating the "head" 
structure. At larger distances they fan out in a char- 
acteristic tail that extends over many tens to hundreds 
of kpc. When the flow of intracluster medium (ICM) 
impacting these galaxies (in their rest frame) is super- 
sonic (Mach number larger than 1), the ram pressure o f 
the ICM causes the jets to bend (jBegelman et al.lll979[ ). 
If the flow is transsonic (Mach number ~ 1), the ther- 
mal pressure gradient of the interstellar medium of these 
galaxies, due to the ir motion through th e ICM, deter- 
mines the bending ([ Jones k, Owenl 11979). In this last 
model, the inflow is decelerated and heated by a bow 
shock in front of the galaxy, which also generates a tur- 
bulent wake that re-accelerates the relativistic particle 
population in the tail and illuminates the tail. 

The radio contours of IC 310 show an extended emis- 
sion, pointing away from the direction of NGC 1275. The 
length of t his tail measured in radio varies between 0.12 ° 
and 0.27° (|Siibring fe de Bruvnl[l99l lLal fe Raoll2005[ ). 
The X-ray image of IC 310 observed by XMM-Newton is 
compatible with a point-like emission from the core and 
with no X-ray emissio n from its ext ended radio struc- 
ture (|Sato et al.ll2005D . Interestingly. iSato et all (|2005f ) 
also showed that the X-ray emission may originate from 
the central AGN of the BL Lac-type object. Other ob- 
served characteristics of IC 310 (e.g. no strong emission 
lines, spectral indexes in radio and X-ray) s uggest that it 
may also be a dim (weakly beamed) blazar (|Rector et al.l 
[1993) . 

T he LAT instrument on board the Fermi satel- 
lite (lAtwood et all 120091 ) has recently detected IC 310 
(jNeronov et al.ll2010D with 5 (3) photons above 30 GeV 
(above 100 GeV). At lower energy (i.e. from 100 MeV to 
1 GeV) , only the cl uster' s central galaxy, NG C 1275, is 
visible (jAbdo et al.ll2009l: iNeronov et al.ll2010l) . 

In this letter, we present the results of recent obser- 
vations of the Perseus cluster performed with the Major 
Atmospheric Gamma Imaging Cherenkov (MAGIC) tele- 
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scopes and the detection of IC 310 (ATel #2510). 

2. OBSERVATION AND ANALYSIS 

MAGIC consists of two 17 m Imaging Atmospheric 
Cherenkov Telescopes (IACT) located at the Roque 
de los Muchachos, Canary Island of La Palma (28°N, 
18° W), at the height of 2200 m a.s.l. The commission- 
ing of the second MAGIC telescope was finished in 2009, 
and since the en d of 2009 both teles copes work together 
in stereo mode (|Cortina et al.ll2009D . The stereo obser- 
vations provide an excellent sensitivity of < 1% of the 
Crab flux (C.U. ) 1 in the medium energy range in 50 hr 
of observations ( Colin e t al. 200^). 

The MAGIC I telescope was used to observe the 
Perseus cluster for a total of 94 hr between Novem- 
ber 2008 and February 2010. The analy sis of 33.4 hr 
of data taken in 2008 was presented in lAleksic et al.l 
(|2010al) ; it focused on the physics of the P erseus cluster 
and NG C 1275. The skymap presented in lAleksic et al.l 
(|2010al) does not show significant excesses in the IC 310 
position 2 . Since the end of October 2009 the second 
MAGIC telescope was also taking data, allowing the 
stereoscopic analysis. 

Obs ervations were per formed in the so-called wobble 
mode (jFomin et al.lll994[) , with data equally split in two 
pointing positions offset by 0.4° from the direction of 
NGC 1275. IC 310 was in the field of view at the angular 
distance of 0.25° and 1° for individual wobble positions. 
Since the 7-ray collection area in the latter case is signif- 
icantly lower (by a factor of ~ 3), only data with 0.25° 
offset were used for the signal search and for obtaining 
the spectrum and the light curve. 

We separately analyzed the MAGIC I single telescope 
(hereafter mono) and stereo data. The mono and stereo 
are only partially independent systems, differing in the 
analysis method, thus there can be some residual sys- 
tematic error between them. For single telescope ob- 
servations, the system atic error in the d etermination 
of the flux is ~ 30% (| Albert et al.ll2008al) . Addition- 
ally, IC 310 was not observed in the standard wob- 
ble observation mode; this increases the systematic er- 
ror from the background estimate. The mono data 
are treated with the standard MAGI C I analysis chain 
(jAlbert et all l2008at I Abu et al.l I2009D . In the analysis 
of the stereo data we took advantage of the impact pa- 
rameters with respect to each telescope and the height 
of the shower maximum. Those stereo parameters im- 
prove the ga mma/hadron separation and the energy re- 
construction (|Aharonian et al]|1997t MAGIC stereo per- 
formance paper, in prep.). We used a new method to re- 
construct the arrival directions that improves the angular 
resolution and the sensitivity of the MAGIC telescopes. 
The method is based on the "DISP RF" technique (see 
lAleksic et al.l l2010bl for details), adapted to the stereo 
observations. The estimate of an arrival direction is per- 
formed independently for each telescope using the shape 
and time information of a corresponding image. In order 

1 Hereafter C.U. stands for Crab Units, defined as the frac- 
tion of the Crab Nebula flux defined by Eq. 2 of lAlbert et al.l 
(2008a), that corresponds, e.g. for energies above 300 GeV, to 
12.4 X 10" 11 cm' 2 s _1 . 

2 Throughout this paper, we refer to the IC 310 position as 
the p osition reporte d by the NASA Extragalactic Database (NED) 
from |2M ASS (200!) catalog 
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Fig. 1. — Fermi-hAT count rate image of the IC 310 region 
in three energy bands: 0.3 — 3 GeV (red), 3 — 30 GeV (green) 
and 30 - 300 GeV (blue). The position of IC 310 and NGC 1275 
are shown with green crosses. Solid (dashed) circles show the sig- 
nal (background) integration regions used for obtaining the Fermi- 
LAT spectrum. 

to obtain the best performance, those are combined with 
the crossing point of the main axes of the images from 
both telescopes. 

We also analyzed Fermi data taken during the period 
between 20 08 August 4 and 201 July 15 following the 
approach of lNeronov et al.1 (]2010l ) . The data were filtered 
with gtselect tool and the Fermi exposure at the source 
position was calculated using gtexposure tool 3 . Only 
events of "diffuse" class were retained in the analysis. 
We obtained the spectrum of IC 310 in the 2 - 200 GeV 
energy band in two different ways. First, we performed 
a spectral fit using the standard Fermi unbinned likeli- 
hood analysis, taking i nto account all so urces from the 
I s * year Fermi catalog (|Abdo et all 120101 ) within 10° of 
IC 310. We then obtain spectral points by extracting 
photon counts from the circle of radius 0.3° centered 
on the source. Taking into account the proximity to a 
bright, nearby source (radio galaxy NGC 1275), we esti- 
mated the background by taking three apertures at 0.6° 
away from NGC 1275 (see FigO). 

3. RESULTS 

After the data quality check, we obtained a sample 
for the period from October 2009 to February 2010, of 
20.6 hr of MAGIC stereo data. The theta 2 (squared dis- 
tan ce between true and reconstructed source position; see 
e.g. iDaum et al.lll997t) distribution of the signal coming 
from IC 310 and the background estimation are shown in 
Fig. [2] (left panel). We found an excess of 106 events, cor- 
responding to a 7.6 a significance (c alculated according 
to the prescription bv lLi fe Malll983l Eq.17). 

The source was also detected in the 27.5 hr of mono 
data (September 2009 - February 2010) with a signifi- 
cance of 8.6 a. Note that since part of the MAGIC I 
data set are also used in the stereo analysis, the two 
significances are not completely independent. The cor- 
responding alpha (the angle between the main axis of 
the 7-ray induced shower image and direction to the 

3 See http://fcrmi.gsfc.nasa.gov/ssc/data/analysis/scitools/ for 
details on the Fermi analysis, we used the software version vyrl5p2. 



true source position; see e.g. I Albert et al.ll2008al) distri- 
bution is also shown Fig. [5] (middle panel). The differ- 
ent signal significance obtained in stereo (7.6er) is sim- 
ilar to the one of mono scaled to the same observa- 
tion time (8.60- x ^20.6 hr/27.5 hr = 7.5<r). This is 
because the mono data have been taken over a longer 
time period, including a higher emission state in Octo- 
ber 2009 (see bel ow) . Moreover, in significance calcu- 
lated according to iLi fe Mai (|1983l ) the background is es- 
timated using both ON and OFF measurements. Thus 
for a high signal/background ratio (as in the case of 
excellent gamma/hadron separation obtained in stereo 
observations), the background is overestimated and this 
lowers the significance. On the other hand a simple cal- 
culation of excess/ a/OFF scaled to the same observation 
time gives a higher value for stereo, 10.2, than for mono, 
9.2, observations. 

It is interesting to note that the 11.2 hr of good qual- 
ity, mono data taken at the end of 2008 do not show 
any significant excess at the position of IC 310 (see 
Fig. [21 right panel). These data yield an upper limit 
for the flux F(> 300GeV ) < 1.9% C.U. (calculated us- 
ing the IRolke et al.|[2005l method with 95% of confidence 
level and assuming 30% systematic error in absolute flux 
level) . 

At energies above 400 GeV the MAGIC telescopes 
working together in stereo have a point spread function 
(PSF) of ~ 0.06°, defined as a 40% containment radius, 
corresponding to a a of a 2-dimcnsional Gaussian. In 
Fig. |31 we show the significance map of the Perseus clus- 
ter region above 400 GeV. The bright spot is consistent 
with the position of the IC 310. In the panels inserted in 
Fig. [31 we also s how archival (non-simultane ous) IC 310 
VLA radio data (jCondon fc Broderickj[l988l ). 

The MAGIC stereo observations reveal a flat spectral 
energy distribution (SED) between 150 GeV and 7 TeV 
without any visible curvature or cut-off (see Fig.@|. The 
differential flux in units of cm~ 2 s _1 TeV _1 is well de- 
scribed (x 2 / n dof = 2.3/4) by a pure power law: 

dN/dE = (1.1 ± 0.2) x 10- 12 (£7TeV)- 2 - 00±al4 . (1) 

The mean 7-ray flux above 300 GeV obtained from the 
stereo observations between October 2009 and February 
2010 is (3.1 ± 0.5) x 10~ 12 cm~ 2 s~\ corresponding to 
(2.5 ± 0.4)% C.U. Comparing this with the upper limit 
from the 2008 data suggests variability of IC 310 on a 
one- year time scale. 

The light curves of IC 310 7-ray emission above 
300 GeV obtained both with the mono and stereo data 
are presented in Fig. [5] Hints of variability can be seen 
in the data. Fitting the individual light curves assuming 
constant flux yields x 2 / n dof = 27.6/7 (for mono, corre- 
sponding to 3.5cr ) and 17.5/4 (for stereo, corresponding 
to 3. Oct). The largest deviations from the mean value are 
for the intervals 13 — 14 October 2009 (3.1 a in mono), 
and 9 — 16 November 2009 (3 a in mono, 3.2 a in stereo). 

Until February 2010, the Fermi-LKY instrument ob- 
served only three photons wi th energies above 10 GeV 
from the direction of IC 310 (jNeronov et al.H2010f ). It is 
interesting to note that one of those 7-rays was observed 
on 15th of October, nearly coincident with the higher 
flux seen in mono. 

The standard Fermi likelihood analysis gave a "Test 
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Fig. 2. — Theta 2 distribution of the IC 310 signal and background estimate from stereo observations taken between October 2009 and 
February 2010 (left panel). Alpha distribution from mono observations taken between September 2009 and February 2010 (middle panel), 
and November-December 2008 (right panel). Only the pointing position 0.25° away from IC 310 is used. The cuts result in an energy 
threshold (defined as a peak of the differential MC energy distribution) of ~ 260 GeV for both mono and stereo data (see the text for 
details). 
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Fig. 3. — Significance skymap from the MAGIC stereo observa- 
tion (42 hr; both pointing positions) for energies above 400 GeV. 
We also show an enlargement of the IC 310 re gion overlaid with the 
NVSS (NRAO VLA Sky Survey at 1.49 GHz; ICondon fc Broderickl 
1198811 contours (top left inserted panel) and the corresponding 
NVSS image (bottom left inserted panel). The NVSS data were 
obtained with Aladin (Bormarcl et al. 2000). Positions of IC 310, 
NGC 1275 and NGC 1265 are marked with black, white, and green 
crosses respectively. 

Statistics" value of 79 from IC 310 above 1 GeV 
(corresponding to a ~ 9 a detection) 4 . Assum- 
ing a simple power law for the spectrum, we ob- 
tain the differential flux dN/dE = (9.5 ± 2.9) x 
10" 9 (£;/10GeV)- 1 - 58±0 - 25 cm- 2 s- 1 TcV- 1 . The Fermi 
spectral index is very hard, mostly due to the last point. 

4. CONCLUSIONS AND DISCUSSION 

The MAGIC telescopes have detected VHE 7-ray emis- 
sion with high statistical significance from the direction 
of IC 310 with a mean 7-ray flux above 300 GeV, between 
October 2009 and February 2010, of (2.5 ± 0.4)% C.U. 
The source seems to be variable on the time scales from 
weeks to about one year. No significant excesses were 
observed in the 2008 mono data and flares are possibly 
present at the > 3 a level in October and November 
2009. The stereo observations yield a spectrum that is 
well fitted by a simple power law with a spectral index 

4 Note that t his significance is not calculated according to the 
ILi fc Mai 1119831 ) method, see Fermi analysis description (foot- 
note [3j . 
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Fig. 4.— SED of IC 310 obtained with 20.6 hr of the MAGIC 
stereo data (full circles). Open triangles show the flux measure- 
ments from the Fermi-LAT from its first two years of operation . 
Archival X-ray {S ato et al. 2005 ). optical jZwickv fc Kowaliri968 |; 
| de Vauc ouleurs et al. 1991), IR (2MASS 2003; Kiiap p"eTainlT98a 
IRAS 1 119881) and radio [[Gregory fc C ondon 1991; Bcck cretaLl 
19911; ICondon et all [2002; White & Bcckci] IT9921; IDouglas et all 
19961 ) data obtained from the NED database are shown with grey 



dots. The solid line shows a power law fit to the MAGIC data, and 
the dotted line is its extrapolation to GeV energies. We also show 
a zoom-in of the MAGIC points. 
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Fig. 5. — Light curve (in 10-day bins) of the 7-ray emission 
above 300 GeV obtained with the mono (open squares) and the 
stereo (full circles) MAGIC data. The open square with an arrow 
is the upper limit on the emission in November-December 2008. 
Vertical grey lines show the arrival times of > 100 GeV photons 
from the Fermi-LAT instrument. The horizontal dashed line is a 
flux level of 2.5% C.U. 
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of -2.00 ±0.14. 

The MAGIC angular resolution is not sufficient to de- 
termine the location of the VHE emission region within 
the radio galaxy. Therefore, it is not clear whether the 
observed 7-ray emission is connected with the tail of the 
source or if it is produced at the base of the jet, close 
to the central engine of the source (as in blazars). The 
strong indications of variability disfavor the 7-ray pro- 
ductio n at the bow shock, discussed by iNeronov et al.l 
(|2010ft . because in this case the emission should be steady 
on time scales of thousands of years. Variability with a 
time scale of a year (a week) constrains the size of emis- 
sion region to be < 10 18 cm (< 2 x 10 16 cm) across, 
assuming no Dopplcr boosting of the flux, which is much 
smaller than the total size of the tail (~ 10 24 cm). Addi- 
tionally, one can estimate the mass of the central black 
hole, Mbh, of an active galaxy using the correlation be- 
tween black hole mass and the central velo city disper- 
sion of the host galaxy (jTremaine et al.ll2002ft . The mea- 
sured velocity dispersion in IC 310 (230 km/s. [McElrovl 
11995ft yields M BH = 2.4 x 10 s M , corresponding to a 
Schwarzschild radius of Rsh — 7 x 10 13 cm. This indi- 
cates that the most probable location of the 7-ray emis- 
sion regio n is in the innermos t part of the jet (as e.g. for 
M 87, see lAcciari et al.ll2009l) . 

The extrapolation of the IC 310 spectrum obtained 
with the MAGIC telescopes is in good agreement with 
the Fermi spectrum below 60 GeV. On the other hand, 
there is a large deviation in the last energy bin measured 
by Fermi (see Fig. [4]) . The 7-ray flux from our obser- 
vations in this energy bin predicts 0.6 photons, while 4 
photons were observed by Fermi-h AT . Assuming a Pois- 
son distribution, the probability of obtaining > 4 pho- 
tons is 3.4 x 10~ 3 (corresponding to 2.7 standard devia- 
tions) so the discrepancy may be a statistical fluctuation. 
However, the Fermi-LAT, having an energy resolution of 
~ 10%, observed 3 of the 4 photons with nearly the same 
energies (98.5, 105, and 111 GeV). If confirmed by future 
observations, these events may indicate the presence of a 
peculiar peak or bump in the IC 310 spectrum given that 
the remaining Fermi-LAT data agree with the MAGIC 
measurement. The detection of such a relatively nar- 
row feature in the spectrum of a radio galaxy may be 
a clue regarding the particle acceleration mechanism at 
the base of AGN jets (e.g. "direc t" gamma-ray emission 
durin g acceleration of particles, INeronov fc Aharonianl 
[20071) ■ This remains a suggestion since the source seems 
to be variable, and the MAGIC and Fermi data used 
here were not taken simultaneously. 

The combined MAGIC and Fermi spectrum (besides 
the above mentioned bump) is consistent with a flat E~ 2 
spectrum stretching without a break over more than 3 
orders of magnitude in energy (2 GeV - 7 TeV). This 
is similar to the flat VHE spectra of M 87, another ra- 



dio galaxy detected at TeV energies ( Aharo 
[27MlAlbert"d~aT] 2008b; Acciari et al.l 12008ft 



Aharonian et al.1 
Such an 

extended, flat E~ 2 spectru m is hard to obtain in a 
simp le one-zone SSC model (|ReesHl967t iMaraschi et~aT1 
11992ft . Instead, a viable model of emission might be in- 
verse Compton scattering of external IR photon back- 
ground pho tons from accretion flow or from the inner 
jet fsee e.g. fNeronov fc Aharonianl 1200 7ft . Alternatively, 
a flat spectrum can be produced in the hadronic m od- 
els (e.g. the proton blazar model, iMannheiml 11993ft . In 
more complicated, multi-zone leptonic models, the GeV- 
TeV emission of a few slightly shifted inverse Compton 
peaks can als o emulate a flat spectrum (e.g . spine-sheath 
layer model, iTavecchio fc G hiscllini 2008 ). 

Finally, using the model by iDominguez "eTaTl (I20l0ft . 
we find that the change in the spectrum due to the ab- 
sorption in the extragalactic background light radiation 
field is within the error of the spectral slope. 

IC 310 and NGC 1265 were the first two head- 
tail g alaxy discovered by astronomers (jRvle fc Windraml 
11968ft . Now, IC 310 is the first head-tail radio galaxy 
detected in the VHE 7-rays by ground based telescopes 
as well as by the Fermi satellite. Additionally, it is also 
the first source discovered above 300 GeV by the MAGIC 
telescopes working together in stereo mode. This detec- 
tion is important also for the VHE study of the Perseus 
cluster of galaxies. It may imply a substantial injection of 
high energy particles in the ICM by a non-central AGN. 
Those particles might have important consequences on 
the possibl e cluster 7-ray emissi on due to cosmic-ray ac- 
celeration (jAleksic et al.ll2010aj) . The VHE IC 310 de- 
tection significantly enriches our knowledge of the 7-ray 
universe. 
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